Photon Radiation and Dilepton Production 
Induced by Rescattering in Strong Interacting Medium 
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Using the opacity expansion technique we investigate the photon radiation and dilepton produc- 
tion induced by multiple rescattering as energetic parton jet passing through the strong interacting 
medium. The real photon radiation and dilepton invariant-mass spectrums of the bremsstrahlung 
contribution from an energetic quark jet are presented. The leading contribution of total energy 
loss by photon emmision in medium of a higher energy quark jet is found to be proportional to 
the jet energy and has a linear dependence on the thickness of the nuclear target. The rescattering 
contribution for the dilepton is important only for small value of the invariant-mass and for the not 
so fast jet. The contribution fraction of the dilepton induced by rescattering in medium is found to 
be nearly a constant when the ratio of jet energy to Debye screening mass E/ is large. 

PACS numbers: ll.80.La, 12.38.Mh, 12.38.Bx, 25.75.Nq 



I. INTRODUCTION 

Hard processes are considered as a good tool to study 
the properties of the quark matter produced in ultra- 
relativistic heavy-ion collision because it can probe the 
early stage of the evolution of dense system, during 
which a quark-gluon plasma (QGP) could exist for a 
short period of time. One important aspect of hard 
processes islet's energy loss or jet quenching [ll, 0, H, 
0, S 113, EI, US, due to gluon radiation induced 
by multiple rescattering as energetic parton jet going 
through the strong interacting medium. Jet quenching 
in high energy nucleus-nucleus collisions is predicted to 
lead to strong suppression of both single- and correlated 
away-side dihadron spectra at large transverse momen- 
tum as compared io p -\- p collisions at the same en- 
ergy [11,112, 13). These phenomenon have been observed 
in central nucleus- nucleus collisions in recent RHIC ex- 
periments [li,|il,[ll,|il[ll,[ll[23|. A simultaneous x^- 
fit to both single and dihadron spectra can be achieved 
within their minima in the same narrow range of en- 
ergy loss parameter for two different measurements at 
RHIC [l^. This fact provides convincing evidence for 
jet quenching description. 

In addition to the gluon radiation, an energetic quark 
jet suffering multiple rescattering in medium should also 
induce real and virtual photon bremsstrahlung, where the 
emitted virtual photon may further decay into a dilepton. 
Since the (real or virtual) photon interacts with the par- 
ticles in the collision region only through the electromag- 
netic coupling, its mean-free path is expected to be quite 
large [2 1|. One expects the photons can pass through the 
collision region without rescattering and then carry the 
information of medium at the time they have been pro- 
duced, which may help us analyze the properties of the 
hot nuclear matter created in relativistic nucleus- nucleus 
collisions and provide complementary test of jet quench- 
ing mechanism due to gluon emission off their same par- 
ent jet[2l,[2l,[23. 



The bremsstrahlung photon production in nucleus- 
nucleus collisions has been previously discussed by many 
researchers. B. G. Zakharov studied [24| the induced pho- 
ton bremsstrahlung from a fast quark produced in ^4 -I- A 
collisions due to multiple scattering in QGP and pre- 
dicted the medium- induced photon emission may en- 
hanced the photon production at high pT \n. A + A colli- 
sions by 30% as compared to p+p collisions, which is not 
consistent with the recent PHENIX measurements [25} . 
By using diagrammatic method in thermal field theory, 
Arnold, Moore and Yaffe (AMY) [111 calculated the pho- 
ton and gluon emission rates of an equilibrated, hot QCD 
plasma to the leading order in both Ue and the QCD cou- 
pling gs{T), where the energy loss of an asymptotic jet in 
QGP is investigated by assuming a very high temperature 
to guarantee gs{T) <C 1. Later S. Turbide, C. Gale, S, 
Jeon and G. Moore (TGJM)[23 utilize AMY formahsm 
to study high px pion and photon production including 
bremsstrahlung process in heavy-ion collisions at RHIC 
and LHC. However, their results overestimate the photon 
production at Au + Au wi th \/ snn = 200 GeV measured 
by PHENIX experiments 

Therefore, it is interesting and beneficial to make an 
complementary calculation of photon production in A+A 
collision with different approaches. In this paper we fol- 
low the framework of opacity expansion developed by 
Gyulass, Levai and Vitev (GLV) in Ref. ^] to study the 
bremsstrahlung photon and dilepton production induced 
by multiple rescattering as the parton jet propagating 
in the strong interacting medium. The GLV formal- 
ism has been extensively applied to study the high px 
phenomenon at RHIC and the phenomenological studies 
based on GLV formalism have described different mea- 
surements at RHIC very well [1, [1, [3| ■ In GLV formalism 
the contributions of radiative energy loss in expanded in 
powers of the opacity, which is defined as the averaged 
scattering number, L/X with L representing the thickness 
of the target and A the gluon mean-free path. And it has 
been shown that the opacity expansion series converge 
very rapidly and the first order contribution is dominant. 
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For the comparison of several jet quenching approaches 
including AMY and GLV formalism, please see the re- 
view in Ref . [§| . 

As a first step we concentrate on the production of pho- 
ton and dilepton emitted off a fast quark jet at the first 
order opacity. In our calculation the QCD coupling gs 
from the scattering point where a quark jet interacts with 
a parton in a QGP system, is put in a potential function 
and absorbed into opacity. Therefore, compared with 
amy's study [2^, our result of photon emission rate at 
first order opacity is also to the leading order in both ag 
and Us- Our conclusions about bremsstrahlung photon 
in strong interacting medium are qualitatively consistent 
with amy's studies. We confirm that the leading con- 
tribution of energy loss by photon emission has a linear 
dependence on the thickness of the nuclear target L[28|, 
which is different from the energy loss by gluon radiation 
where Ai5 quadratically depend on the thickness L. We 
observe that the rate of photon radiation as a function of 
energy fraction x carried by radiated photon has a peak 
at small x ~ 0.3, whereas in Zakharov's study the peak of 
the rate by induced photon correction is at a: ~ 0.9. This 
may explain why a larger enhancement factor of photon 
production at high pT in ^ -f ^ is found in Zakharov's 
study. For a jet created in A + A collisions, multiple 
scattering induces photon radiation as well as gluon ra- 
diation (jet quenching). Compared with the case of p-f p 
collisions, induced photon radiation enhances the pho- 
ton yield while jet quenching decreases the pho ton yield 
from jet fragmentation vn A + A collisions]!^, [IS, US]- 
If the former effect is overestimated larger than the lat- 
ter effect, the nuclear modification factor for medium and 
large photon production Raa > 1 will appear in ^-l-A 
collisions, not consistent with the recent PHENIX mea- 
surements (25^ . 

For dilepton production it is shown the rescattering 
contribution is important only for small value of the 
invariant-mass and also for the not so fast jet, and the 
contribution fraction of dilepton induced by multiple 
scattering to approach a constant when the quark jet 
energy E/ ^j, > 30, where E denotes the energy of the jet, 
and jjL the debye screening mass. 

This paper is organized as follows. In Sec. II we 
present the calculation of the jet energy loss by pho- 
ton emission in hot medium and the resulting photon 
spectrum due to real photon radiation. In Sec. Ill we 
investigate the invariant-mass spectrum in the dilepton 
production due to virtual photon radiation. Correspond- 
ing numerical analysis are given in these two sections, 
respectively. In Sec. IV we make a brief conclusion. 
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where x is the energy fraction carried by photon from the 
jet parton. The initial and final four momentum of the 

jet is assumed as 



Pi = [P'^,P ,0_l]. 



Pf 



{l-x)p-^ 



(1 - .T)p^ 



-,qj 



(3) 
(4) 



where q_L and m denote the transverse momentum trans- 
fer and thermal mass of the jet parton in the hot medium, 
respectively. 

As proposed by Gyulassy and Wangpjl, the interaction 
between the jet and the target parton can be modeled by 
the static color-screened Yukawa potential, 

Vn = 27r5(g")t;(q„)e~''i"-''"raJj)®T,„(n), (5) 

w(q 2,2 ' (6) 

where as is the strong coupling constant, p the Debye 
thermal mass of gluon in the hot medium, q„ the momen- 
tum transfer from a target parton n at the coordinates 
Yn — {zmYLn)- (j) and Ta^ {u) are the color matrices 
for the jet and target parton, respectively. As assumed 
in Ref.jJ, all target partons are in the same dr dimen- 
sional representation with Casimir C2(T) {TrTa{n) = 
and Tr{Ta{i)n{j)) = 5^j5abC2{i)dT / (Ia). For color ma- 
trices generators of the jet in the dimensional repre- 
sentation, Tr{Ta{j)Ta{j)) = CndR. 

As shown in Fig.l, we denote A^o, -^1 and A^2 as the 
scattering amplitudes for self-quenching, single rescatter- 
ing and the double Born rescattering, respectively. The 
squared total amplitude is given by 



\Mq + Mi 
|A/(o|' + |Xi|' 
+2Re{M2Ml) 



Mo 



2Re{MiM*o) 



(7) 



For the self-quenching shown in Fig.l we get 

Mo = iJ{p + k)e'^P+'''>-y°'R.o , (8) 
where the radiation amplitude TZq can be expressed as 



II. INDUCED REAL PHOTON RADIATION 
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(9) 



Consider a source localized at yo — (^Oiyo) that pro- 
duces a parton jet described by a wave packet J{p). The 
parton jet undergoes multiple rescattering by exchanging 
a gluon with a target parton and emits a photon with 



g is the coupling constant between quark and photon. 

For multiple rescattering we assume the packet J(p) 
vary slowly over the range of momentum transfers sup- 
plied by the potential, the rescattering centers are well 
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FIG. 1: Feynman diagram for photon radiation induced by 
self-quenching, single and double Born rescattering in the 
medium. 



separated and the distance between two successive scat- 
tering is large compared to the interaction range, 



(10) 



For single scattering illustrated in Fig.l the rescatter- 
ing amplitude of the first Feynman diagram can be ex- 
pressed as 



Ml = iJ{p + k)e 
dq. 



N 

i{p+k)-yo 



2g{l - x)p+e-p^ 



-yoi 



27r [(p-g)2+ie][(p + fc-g)2 



Similarly we get the radiation amplitude for second Feyn- 
man diagram of the single rescattering in Fig.l as 



-i2g{l - x) 



(ki 



a;qi 



N 



iu)o(zj-za) 



(14) 



where N is the number of the targets in the medium. The 
integration over Qz is performed by closing the contour 
below the real axis in the complex qz plane. Two poles 
qz' = —ie and qz" — —LUQ — ie contribute to the scattering 
amplitude. Here loq can be written as 



UJQ = 



+ x^m"^ 



xp'' 



(12) 



The contribution from the pole corresponding to the po- 
tential singularities can be neglected because of the sup- 
pression factor exp[— ^(z — zq)]. From Eg. pTj) we deduce 
the radiation amplitude for first Feynman diagram of the 
single rescattering in Fig.l as 



e± ■ kj 



N 



k^ + x'^m? ^ 



(13) 



Then the total radiation amplitude for single scattering 
is given by 

7^l = Tl\^Tl\=-i2g{\-x) 

N 

X • ^ [b - (B - C) e*"°(^^-^«)] , (15) 



where 
B = 



C = 



xq± 



(k_L — xq±)'^ + x'^m? 



(16) 



For the double Born virtual interaction shown in Fig.l, 
yi = 1/2 = Vj, the rescattering amplitude M2 reads 



X g-i{^i±+^2±)-(y ±j-y ±0) — x)^p'^'^ 



dq2z 
2tt 



X v{q2)Ta,Ta,ij)A{p - q2) J ^« (qi )T,, T^, ( j) 

X e-'(«i^+«^^)("^-"")2g[e ■ (p - qi - 92)] 

X A{p-qi-q2)A{p + k-qi-q2). (17) 

Integrate over qiz by closing the contour below the real 
axis in the complex qiz plane, the residue of the poles 
qiz = —q2z — *e and qiz ~ —q2z — cjJq — ie contribute to 
the integral. The contribution from the potential singu- 
larities can be ignored because Zj — Zq ^ 1/ H- so A^2 
can be rewritten as 



(11) = iJ{p + k)e 



N 

i{p+k)-yo 



g-i(qii+q2i)-(yj_i-yio)^2 — .t)^p+ 
-2ig 



(27r)2 J {2nY 
2 f dq2z 



2-K 



wo(l - x)p^ 



rw(g2z, q2x)w(-'?22, qi±) A(p - 92) 



X e-pi 



1 



^ii^a{zj-za) 



Ta,Ta,{j)Ta,Ta,U) . (18) 



For q2z integration, there is no exponentially suppressed 
factor, so the potential singularities q2z — — i/^i and —ifi2 
contribute to the integral. By closing the integration 
contour below the real axis in the complex q2z plane and 
taking the residue in the poles q2z = — and — 
—i\/qL2± + the radiation amplitude 



N 



n'^{Born) = -5(1 - x)e± • B ^ [l 



^iLL>a{zj ~za) 



(19) 
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In a similar way, from the last Feynman diagram in 
Fig. 1 we obtain the rescattering amplitude Ai2- The 
corresponding radiation is given by 



N 



nl{Born) = ~g{\ - x)€± ■ B ^ 



„iwo(2j-zo) 



(20) 



The total radiation amplitude from the double Born 
virtual interaction can be expressed as 



7^f = 7^2 (Sorn) +7^2 (Born) 
= -g{l-x)Ne^-B. 



(21) 



As shown in Ref. 4] , in opacity expansion technique we 
make following simplification: define the relative trans- 
verse coordinate h±j = y±j — y±o, it varies over a large 
transverse area A± . Then the ensemble average over the 
rescattering center location reduces to an impact param- 
eter average, 



= J d^hj^dzjdzop{zj,zo) . 



(22) 



Along the longitudinal direction, the distribution den- 
sity for the locations of the rescattering center is defined 
by 



9{L-z) 

L/2 
Ojz-zo) 

L/2 



Exp 



Exp — 



L/2 
L/2 



(23) 



Taking the initial color average and the ensemble av- 
erage we obtain 



= / dzodzp{zo,z) 



1 CbC2{T)\v{0,ci^)\^ 



(27r)2 



7^l|^24) 



The elastic cross section with the small transverse 
momentum transfer between the jet and target partons 



CrC2{T) |t;(0,q^)|- 
dA (27r)2 



(25) 



One can defining |u(0, q_L)p as the normalized dis- 
tribution of momentum transfers from the scattering 
centers j^]. 



1 dfj, 



el 



l«(0,q±)P = - 



IT (qi -I- /i2 



where 



(26) 



(27) 



to insure J'^-^'"""' (i2qj^|u(0, q^)p = 1. In numerical esti- 
mation, we take that q^^^^. — oo, ^eff ~ M- 

Using the normalized distribution |ti(0, qj.)^, we can 
rewrite Ea. ([M|) as 



(l-^iP) = f /rfV|S(0,q^)p 

X / c^zoc^zp(zo,^)|7^l|^ (28) 



where 



L Na, 



A A, 



(29) 



n represents the mean number of the rescattering which 
is defined as the so-called opacity Q. 

The first interference term RbIMiMq) in Eq.([7|) 
doesn't contribute to photon radiation probability be- 
cause of Tr{Ta{j)) = 0. The contribution from the dou- 
ble Born virtual interaction is given by 

{2Re{M^M*o)) = jj d\^\v{0,<i^)f 

X / c^zo^^z/9(^o,z)2i?e(7^f7^o).(30) 



Denote N as the number of the radiated photon, we 
have 



dN = {\M\^ 



d^k 



(2nf2uj 



(31) 



The transverse momentum spectrum of the radiated pho- 
ton can be expressed as 



dN 



1 1 



d\k^\^dx 4(27r)2 



{\M\' 



(32) 



To the first order of opacity expansion, the square of 
the rescattering amplitude is given by 

= {\Mo\^ + \Mi\^ + 2Re{MiM*o)) 
L f /i^d^q^ 



|7^ol 



dzodz 



A J n{q]_ +A*^)2 

X p(zo,^)[|7^lp + 2i^e(7^f7^s)], (33) 



where 



|7^lp + 2i?e(7^f7^S) =4.9^(1 -x)2[(B-C)2 

- 2cos(cjo(^-^o))(B2-B-C)]. (34) 

The term with cos{luo{z — zq)) reflects the destructive 
interference arising from the Abelian LPM effect [295 . The 
integration of this term over target distribution gives 



I{lU0,L) = I dzodzp{zQ,z)cOs{LUo{z - Zq)) 

4 



4 + uj^L^ 



(35) 
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Denote dA^o+i = dNo + dNi, Nq is the number of radi- 
ated photon at zero-order opacity expansion correspond- 
ing to self- quenching term |7WoP, Ni the number of radi- 
ated photon at first-order opacity expansion correspond- 
ing to the sum of single rescattering term and 
the interference term between self-quenching and double 
born virtual interaction 2i?e(A^^A^Q). Defining 



u = |q_L 
we obtain 



dNo 4ae (1 - x)^ 



y 



dxdy Qtt x (j/ + x'^wY 



m^lii^ , (36) 



(37) 



dNx 4ae (1 - xf L 



du- 



1 



Tf{x,u,y), (38) 



dxdy Qtt x XJ (l + u)^' 
where function f{x, u, y) is given by 

{y + x'^u){y + x'^w + x'^u) — Ax'^uy 



f{x,u,y) 



^ {{y + x^u + x'^w)'^ — Ax'^uy)'^ 

, ^2 + (1 - ^(^0' ^)) — ^ 

(y -|- x'^wy y + x'^w 



y — X w 



y — x'^w — x'^u 



y + x^w ^[y + x'^u + x'^wY - Ax'^uy 



(39) 



In getting above expressions the integral of the angle be- 
tween k_L and qj_ has been integrated out and we have 
chosen the coupling constant g = 2e/3 for u light quark 
jet. 

If setting the mass of quark to vanish, from Eq. p6p we 
see w = 0, then function f{x, u, y) in Eg. pop reduce to 



f{x,u,y) 



1 



\y — x^u\ 



(l-/(-o,2.))i(l-f^ 
y \ \y- x^u\ 



.(40) 



This function have two coUinear divergences. One is from 
the initial-state bremsstrahlung as y ~ (or kj^ ^ 0) 
along the initial jet direction; the another one is from 
the final-state bremsstrahlung along the final jet direc- 
tion as |kj^| ~ xjq^l (or emission angle 9e is near to the 
scattering angle 0s), here 9e and 9s approximate 



(41) 



k+ 



Above two collinear divergences is caused by setting 
quark jet's mass to be zero. As 9e > 9s, the third term 
in Eq. (j40p vanish, so there is no interference. In purely 
QED case, as 9e < 9s the radiative ener gy l oss of the 
charged particle has been studied in Ref.[30] in which 
the charged particle undergoes multiple rescattering by 
exchanging photon. Here we investigate the transverse 
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FIG. 2: (Color online) The production rates of the photon 
emitted off two quark jets with energy E/fi = 10, 40, respec- 
tively. The upper two curves (the dot for the jet with energy 
E/fx = 40, the solid for the jet with energy E/n = 10) are 
for the contributions to the first order opacity (denoted as 
"0-1-1"). The long-dashed curve is for the self-quenching con- 
tribution (denoted as "0") from the jet with E/fi = 10 while 
the short-dashed curve for the induced contribution (denoted 
as "1") from the same jet at the first order opacity. The 
opacity is chosen as L/A = 6. 
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FIG. 3: (Color online) The radiation rate of the electromag- 
netic energy from a energy-given quark jet by radiating pho- 
tons carrying energy fraction x of the jet. 



momentum spectrum of radiated photon and the corre- 
sponding energy loss of energetic quark jet which under- 
goes multiple rescattering by exchanging gluon. 

To avoid the collinear divergences, in the following we 
take the thermal mass of quark and gluon from "hard 
thermal loops" [31] as 



9'sCfT' 



Nf 



(42) 
(43) 



where Nc and Nf are the number of the color and flavor 
of quark, respectively. After taking into account the ther- 
mal mass of quark, the collinear divergences in Eq. (j39p 
disappear. In following numerical calculation, we will 
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FIG. 4: (Color online) The ratio between the rescattering con- 
tribution at the first order opacity and that of corresponding 
BH-limit for a quark jet with different energies. 

choose fi = 0.5GeV given by Eq. (|43|) . the mean free path 
A — Ifm, and consider the kinetic limits of the photon's 
transverse momentum (lo| . 

< ki < 4:E^x{l - x) . (44) 

Shown in Fig[2] is the production rates of the pho- 
ton emitted off two energy-given quark jets, respectively. 
The total photon production rate ("0-1-1" order opac- 
ity) decrease with the increasing photon energy {xE) 
and increases with the increasing jet energy. The varia- 
tion tcndence of the total rate is the same as AMY and 
TGJM's studies on brcmsstrahlung contribution 26, 27]. 
Shown in Fig[3] is the radiation rate of the electromag- 
netic energy from a energy-given quark jet by radiat- 
ing photons carrying energy fraction x of the jet. The 
variation tendence of the radiation rate is also similar 
to that of amy's studies on brcmsstrahlung contribu- 
tion. Of interest is that our "0" and "1" order results 
can be directly compared with the "uaccMm" and the 
"mdwced" contributions in Zakharov's studies, respec- 
tively. As /i = 0.5GeV is chosen in Eq. (P5| . the jet en- 
ergy shown in Fig. [3] has a value E = 20GeV same as an 
example in Zakharov's studies ^2^. The variation ten- 
dence of our "0" order contribution is similar to that of 
"vaccum^^ contribution in Zakharov's studies. However, 
the peak of "1" order correction contribution is at a; = 0.3 
in our study while the peak of "mduced" correction con- 
tribution is at a; = 0.9 for finite kinematic boundaries 
in Zakharov's studies. For a jet created in A + A col- 
lisions, multiple scattering induces photon radiation as 
well as gluon radiation (jet quenching). Compared with 
the case oi p + p collisions, induced photon radiation en- 
hances the photon yield while jet quenching decreases 
the photon yield from jet fragmentation in A + A colli- 
sions. If the former effect is estimated smaller than the 
latter effect, the nuclear modification factor for medium 
Pt photon production Raa < 1 could be exp ected in 
A-l-A collisions, close to experiment findingsf25j. Our re- 
sult for induced photon contribution smaller than that in 



FIG. 5: (Color online) The total electromagnetic energy loss 
as a function of the energy of a quark jet. 




L/A 

FIG. 6: (Color online) The total electromagnetic energy loss 
of a quark jet {E/ jj, — 5 ~ 100) as a function of the thickness 
of the strong interaction medium. 

Zakharov's study may explain why a too larger enhance- 
ment factor of photon production at high px in A + A is 
found in Zakharov's study. 

It is the same as pointed out in Ref.[l], from the rel- 
ative phase factor ujq(z — zq) in Eq. (j34p we define the 
formation time of the photon radiation as 



When z — Zq ^ T, the photon radiation reach the Bethe- 
Heitler (BH) limit in which the intensity of the induced 
radiation is additive in the number of rescattering. The 
BH-limit corresponds to I{llJo,L) = in Eq. ((35|) . Shown 
in FigU] are the ratio between the rescattering contribu- 
tion at the first order opacity and that of correspond- 
ing BH-limit for a quark jet with different energies, 
E/n = 5 — 100. From FigH] we see that, at fixed jet 
energy, as x approaches small value, r becomes small, 
the BH-limit is approached; at fixed E/fi becomes 
small, T becomes also small (in Eg. (j45)) . p+ ~ "^E), the 
BH-limit is achieved again. 

Integrating over x in FigO with different value of the 



quark jet energy, 

/\E/E ^ j xdx j dk^ 



dN 
rfk^ dx 



(46) 



we obtain the total electromagnetic energy loss as a func- 
tion of the energy of a quark jet, shown in FigO For a 
fixed target thickness L/\ — 6, the self- quenching contri- 
bution increases with the increasing jet energy while the 
rescattering contribution ("1" order opacity correction) 
is almost a constant /S.E/E/ae — 0.1 when E/ ^ > 5GeV. 
It show that for a high px quark jet the leading contri- 
bution of the total induced electromagnetic energy loss is 
proportional to the jet energy. Shown in Fig[6]is the total 
electromagnetic energy loss of a quark jet {E/fi = 5—100) 
as a function of the thickness of the strong interaction 
medium. It is clear that the leading contribution of en- 
ergy loss by photon emission has a linear dependence on 
the thickness of the nuclear target, which is different from 
the energy loss by gluon radiation where AE quadrati- 
cally depend on the thickness [l^]. Similar results 
are also observed in Ref . I28l . 






III. THE INDUCED DILEPTON PRODUCTION 

What we study in above section is the real photon ra- 
diation when an energetic parton jet produced in heavy 
ion collisions propagates inside the strong interacting 
medium. If the radiative photon induced by multiple 
rescattering in above process is virtual, the virtual pho- 
ton can decay into a lepton pair (U). In this section we 
investigate the induced dilepton production. 

The Feynman diagram for dilepton production induced 
by self-quenching, single and double Born rescattering in 
the medium is shown in Fig[6l The amplitude for the 
dilepton radiation reads 



M 



C, 



fc2 



-u{l'^){~ie^i,)u{l ), 



where 



(47) 
(48) 



(49) 
(50) 



and is relate to the radiation amplitude for the virtual 
photon. 

The invariant-mass spectrum for dilepton production 
can be expressed as 



dN 



dAP SnlVP 



(27r)32fcO ' 



(51) 



where = , M is the invariant mass of the dilepton. 

For virtual massive photon field, we can choose its po- 
larization vector e to satisfy e ■ k = 0. Correspondingly 



FIG. 7: Feynman diagram for dilepton production induced 
by self-quenching, single and double Born rescattering in the 
medium. 



we have [321 



^e^(A)e*(A) = -5^^ 



(52) 



Then we obtain 



2_ Mi^f 

IVP 



(53) 



In getting last equality the Ward Identity has been used. 

As an approximation, we consider a massless quark jet 
and assume p+ fc"*" 3> M . From Eqs.fSTjl and (fSS]) we 



see that e^C^ is the radiation amplitude for a on-shell 
virtual photon with mass M reduced by rescattering. Its 
calculation is very similar to the calculation of the ra- 
diation amplitude for a massless real photon in Sec. II, 
the only difference is that the light-cone 4-momentum in 
Eq.([T]) should be replaced by 



KP 



xp 



xp^ 



-,kj 



(54) 



For self-quenching, single and double Born rescattering, 
the radiation amplitude of the virtual photon with mass 
M can be deduced as 



TZ'o = 2{l-x)ge± B' 



(55) 
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7^'l = -2ig{l-x) 

N 



where 
B' 

and 



+ M2 



c 



xqj_ 



M2 



UJn 



xp~' 



(56) 



(57) 



(58) 



(59) 



Correspondingly, the square of the radiation amplitude 
for the virtual photon to the first order in opacity expan- 
sion can be expressed as 



qi+/z2)2 



dzodzp{zo, z)[|7^'lp + 2i?e(7e'f 7^'o)] (60) 



where 



2i?e(7^f O = V(l - x)2[(B' - C')2 
-2 cos(tj^,(z - zo))(B'^ - B' • C')] (61) 



It is the same as in Eg. (p4|) . the term with cos{uj'q{z — zq)) 
represents the destructive interference resulting from the 
Abehan LPM effect 29]. 

Substituting Eq. (f60|) into Eq. (f5T|) . we can obtain the 
invariant-mass spectrum arising from self-quenching at 
zero order in opacity, 



4a2 {l-xy 



dNp _ 
dxdydw 277r2 xw {y + wY 



(62) 



and the invariant-mass spectrum arising from single and 
double Born rescattering at first order in opacity. 



dNi 
dxdydw 

where function 
f{x,u,y,w) 



Aal {\-xfL 
A 

^f{x,u,y,w), (63) 



277r2 xw X 
1 



du- 



(l + u) 

{y + x'^u){y + w + x'^u) — Ax'^uy 
^ {{y + x^u + w)2 — 4x^uy)^ 



(y + w)2 
( y-w 



y + w 
y — w — x'^u 



y + w ^[y + x'^u + wY - Ax'^uy 



(64) 
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FIG. 8: (Color online) The invariant-mass spectra ("0-1-1" 
order) of the dilepton emitted off quark jets with different 
energies, E/ p = 5, 10, 20, 40, 60, 80, 100, passing through the 
medium with the thickness, L/A = 6. 




FIG. 9; (Color online) The comparison of the induced con- 
tribution (denoted as "1 order") and the self-quenching con- 
tribution (denoted as "0") of the invariant-mass spectra of 
the dilepton emitted off an energy-given quark jet, E/ p — 40, 
passing through the medium with the thickness, L/X = 6. 



. 8 



L/A=6 

M/^=0 . 001,0 .01,0.1,1 




FIG. 10: (Color online) The ratio of the invariant-mass spec- 
tra between 1-order and (0+l)-order in opacity expansion. 
The invariant-mass is chosen as M//i = 0.001,0.01,0.1,l, re- 
spectively. 
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and 

^-Iq^IVm', y = \h^\V^i\ w^M^/t?. (65) 

Shown in FiglS] is the total invariant-mass spec- 
tra ("0-1-1" order opacity) of the dilepton emit- 
ted off quark jets with different energies, E/^m = 
5, 10, 20, 40, 60, 80, 100, passing through a target with the 
thickness, L/X = 6. The invariant-mass spectrum de- 
creases with the increasing invariant-mass while increases 
with the increasing jet energy. Shown in Figl9] is the 
comparison of the induced contribution ("1 order" opac- 
ity) and the self-quenching contribution ( "0" opacity) of 
the invariant-mass spectra of the dilepton emitted off an 
energy-given quark jet. For a jet with energy E/ ^ = 40, 
the rescattering contribution is very small and can be 
neglected in the large invariant-mass region (especially 
when M//i > 1). In order to know clearly how much 
the rescattering contribution is in the total contribution, 
we plot FigHni to check the ratio of the invariant-mass 
spectra as a function of the jet energy between 1-order 
and (O-l-l)-order in opacity expansion. The invariant- 
mass is chosen as M//i = 0.001,0.01,0.1, 1, respectively. 
The plot demonstrates that the dilepton production in- 
duced by rescattering is important for small value of the 
invariant-mass and also the jet energy. In addition, the 
contribution fraction by rescattering is found to be nearly 
a constant when E/ ^ > 30. 

IV. CONCLUSION 

In this paper we apply GLV's opacity expansion tech- 
nique to investigate the induced photon radiation and 
dilepton production in strong interacting medium. The 
real photon radiation and the dilepton invariant-mass 
spectrums are presented to the first order in opacity 
expansion. Our study focuses on the electromagnetic 
bremsstrahlung of an energy-given quark jet passing a 
thickness-given QGP system and the analytical deriva- 
tion is combined with numerical calculation. 

We show that for real photon radiation due to Abelian 
LPM effect the parton jet's energy loss has a linear 



dependence on the thickness of the targets instead of 
quadratic dependence arising from non- Abelian LPM ef- 
fect for gluon radiation. And the leading contribution 
of the total induced electromagnetic energy loss of a fast 
quark is proportional to the jet energy. 

The dilepton production by the decay from the vir- 
tual photon off a quark jet induced by multiple rescat- 
tering in medium has been investigated. It is shown that 
the invariant-mass spectrum decreases with the increas- 
ing invariant-mass while increases with the increasing jet 
energy. The rescattering contribution is important for 
small value of the invariant-mass and also for the not so 
fast jet. In addition, the contribution fraction by rescat- 
tering is found to be nearly a constant when the quark 
jet energy E/^ > 30. 

Of course, the study presented here is just the first 
step for photon and dilepton production in medium. To 
make a complete study of photon and dilepton produc- 
tion in A + A collisions many other contributions should 
be taken into account, such as the prompt photon pro- 
duced by the initial hard scattering, fragmentatedpho- 
tons by jet fragmentation, jet-photon conversion [33j, et 
al. A systematical study will also include other nuclear 
effects such as Cronin effect, ( ant i-) shadowing effect and 
final jet quenching effect for a jet fragmentated into a 
photon, just like our current studies [H, [11] for high pT 
hadrons in A + A collisions. This kind of detailed cal- 
culating and analyzing for large transverse momentum 
photon production will be left for a future's study. 
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